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Abstract 

Pomeron loop effects (gluon number fluctuations) on deep inelastic scattering (DIS) in the 
fixed coupling case are studied. We show that the description of the HERA data for the inclusive 
structure function F2(x,Q 2 ) for x < 10 -2 and 0.045 < Q 2 < 45GeV 2 is improved once pomeron 
loop effects are included. 
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1. Introduction 

Recently, there has been a tremendous theoretical progress in understanding the high 
energy QCD evolution beyond the mean field approximation, i.e beyond the BK-equation [1] 
It has been understood how to include the fluctuations of gluon numbers (Pomeron loops) 
in small-x evolution [2-4]. A relation between high-energy QCD and statistical physics 
was found which has clarified the physical picture and the way to deal with the dynamics 
beyond BK-equation. After including fluctuations, the evolution becomes stochastic and 
the physical amplitude, which one obtains by averaging over all individual realizations. 

2. Mean field approximation 



The Balitsky-Kovchegov (BK) equation [1] gives the evolution with rapidity Y = 
ln(l/a;) of the T-matrix for a qq dipole scattering off a target which may be another 
dipole, a hadron or a nucleus[5], 

?L~a s [T-TT} (1) 

One of the hallmarks of the BK-equation is the so called geometric scaling behavior of 
the scattering amplitude, 
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T(r,Y)=T(r'Qi(Y)) 



(2) 



where r = \x± — y±\ is the size of dipole and Q S {Y) is the saturation momentum defined 
such that T(r = 1/Q S ,Y) be of order 1. The T-matrix shows the geometric scaling 
behavior, i.e., it depends only on a single quantity r 2 Q 2 s {Y) instead of depending on r 
and Y separately which seems well supported by HERA data[6,7]. 

Another important result extracted from the BK-equation is the rapidity dependence 
of the saturation scale, 



Ql(Y) 



2 exp[ax'{X s )Y] 



[aY] 



2(1-A ) 



(3) 



where x(A) = 2^>(1) — -0(A) — ip(l — A) is the eigenvalue of the BFKL kernel. 
Iancu, Itakura and Munier [7] have used the following ansatz for the T-matrix, 



T UM (r,Y) = { 



N 



1 - exp [-aln 2 (brQ s (x))] , rQ s (x)>2 

ln(2/rQ 3 (x)) \ 

1 , rQ s (x)<2, 



(4) 



which obviously includes the features of the solution to the BK equation, to compare 
the theory in the mean field approximation with the DIS data. They have used for the 
saturation momentum the leading Independence of Eq(3), Q s (x) — (x n /x) x , however, 
with A and x being fixed by fitting the DIS data. The constant n — x"(A s )/x'(A s ) w 9.9 
is a LO result coming from the BK-equation, Ao is a constant around 0.5 and a and b 
are determined by matching the two pieces in Eq.(4) at rQ s = 2. 
The "BK-diffusion term" in the IIM-ansatz (4), 



^^^^ 



= exp 



ln 2 (4/r 2 Q2 (x)) 



2 k XY 



(5) 



does explicitly violate the geometric scaling behaviour. We wish to emphasize here that, 
as also shown in [7] , this violation seems required in order to get an accurate description 
of the DIS data. Without it, even allowing A s to be an additional fitting parameter, one 
can not get a better description of the DIS data. For further details on the importance 
of the diffusion term see Ref. [7] . 



3. Beyond the mean field approximation 

3.1. Event-by-event scattering amplitude 



In this paper we will use the Golec-Biernat, Wiisthoff (GBW) [8] and Iancu, Itakura, 
Munier (IIM) [7] models for the event- by- event amplitude. 
The GBW model 



T GBW (r,x) = l-exp 



r 2 Ql{x) 



(6) 
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model /parameters 


x 2 


X 2 /d.o.f 


x (xitr 4 ) 


A 


fl(fm) 


D 


T GBW 


266.22 


1.74 


4.11 


0.285 


0.594 





( T GBW^ 


173.39 


1.14 


0.0546 


0.225 


0.712 


0.397 



Table 1 

GBW model: The parameters of the event-by-event and of the physical amplitude. 

is one of the most simple models which shows geometric scaling, T(r,x) — T{r 2 Q 2 {x)), 
and leads to a quite successful description of the HERA data. 



model / parameters 


x 2 


X 2 /cl-°.f 


x (xl0~ 4 ) 


A 


fl(fm) 


D 


T IIM 


150.45 


0.983 


0.5379 


0.252 


0.709 





^IIM ^ 


122.62 


0.807 


0.0095 


0.198 


0.812 


0.325 



Table 2 

IIM model: The parameters of the event-by-event and of the physical amplitude. 



3.2. Physical scattering amplitude 

Based on the relation between high-energy QCD evolution and reaction-diffusion pro- 
cesses in statistical physics [4,9], the physical amplitude, T(r,Y), is then given by aver- 
aging over all possible gluon number realizations/events 

(T((p - p,(Y))) = J d Ps T(p - Ps {Y)) P( Ps (Y) - (p s (Y))) , (7) 

where have used p = ln(l/r 2 Q§) and p s (Y) = ln(Q 2 s (Y)/Ql), T{p- p s {Y)) is the ampli- 
tude for the dipole scattering off a particular realization of the evolved proton at Y and 
the distribution of p s (Y) is a Gaussian [10]: 



P (Ps) ^ ^=CX P 



(Ps - {Ps)f 



(8) 



The expectation value of the front position, (p s (Y)), increases with rapidity as (p s (Y)) = 
\yl(Q 2 s {Y) / Qq) at high energy [4]. The dispersion of the front at high energy increases 
linearly with rapidity, 

a 2 = 2[(pl)-(p s ) 2 ] =DY (9) 

where D is the diffusion coefficient, whose value is known only for a — > (asymptotic 
energy) [9,2]. Since the values of D and the exponent A of the saturation scale, Q 2 (x) = 
1 GeV 2 (x /x) x , arc not known for finite energies, in what follows we will treat them as 
free parameters. 

4. Numerical results 



We use the dipole picture [11], which is a factorization scheme for DIS, to study the 
proton structure function, which leads to the following expressions for the F 2 structure 
function: 



3 



F 2 (x, Q 2 ) = — f (a T (x, Q 2 ) + a L (x, Q 2 )) (10) 

47T Ct em 

ctt,l(x,Q 2 ) = J dzd 2 r \ip T ,L(z,r,Q 2 )\ 2 a d i P (x,r) (11) 

We fit the ZEUS data [12] for the F 2 in the kinematical range x < 1CT 2 and 0.045 GeV 2 < 
Q 2 < 50 GeV 2 . The interesting ingredient for us in Eq. (11) is the dipole-proton cross 
section, Gdip — 2ttR 2 (T(r,x)), with 2ttR 2 being the outcome of the integration over the 
impact parameter. 

Now, using the GBW model as an event-by-event amplitude, we include the effect 
of gluon number fluctuations by averaging over all events via Eq. (7). The resulting 
(T GBW (r, x)), which breaks the geometric scaling, leads to a better description of the 
HERA data, as can be seen from the comparison of the x 2 values in Table 1. The large 
improvement after including fluctuations seems to indicate that violations of geomet- 
ric scaling, and probably even gluon number fluctuations, are implicated in the HERA 
data [13]. 

In the IIM case fluctuations do improve the description of the HERA data, however 
not much, as can be seen from the comparable x 2 /d.o.f values for T IIM and (T IIM ) in 
Table 2. This is so because the IIM model does already contain the geometric scaling 
violations via the BK-diffusion term, ln(4/r 2 Q 2 ^/v / 2KAF, and describes accurately the 
HERA data, before including fluctuations. Since in the case of the IIM model the fluc- 
tuations do not improve much the description of the HERA data, one may conclude 
that the BK-equation alone should describe the HERA data and that fluctuations may 
be negligible in the energy range of the HERA data. The intention of this work was to 
illustrate the possibility that fluctuations may be present in the HERA data [13]. 
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